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O ■ Abstract 

r^ • We study the strangeness production in antiproton-nucleus collisions at the beam momenta 

Ph ■ from 200 MeV/c to 15 GeV/c and in p annihilation at rest within the Giessen Botzmann-Uehling- 



Uhlenbeck (GiBUU) transport model. The GiBUU model contains a very detailed description of 
underlying antinucleon-nucleon cross sections, in-particular, of the strangeness production chan- 



fj ■ nels. We compare our calculations with experimental data on A, K^ and A production in pA. 

P5 . collisions and with earlier intranuclear cascade calculations. The contributions of various partial 



channels to the hyperon production are reported and systematic differences with the experiment 
are discussed. The possible formation of bound A- and AA-nucleus systems is also considered. 



(N 

> 

^f^ . Finally, results on the H hyperon production are presented in view of forthcoming experiments 



with antiproton beams at FAIR. 
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I. INTRODUCTION 



The production of strange particles in antiproton annihilation on nuclei has been a chal- 
lenge to theory since about two decades when the A and K^ yields and spectra were measured 



in the series of bubble chamber experiments 



3|. The most intriguing experimental fact 



is that the yield ratio A/K^ is very large (> 2) not only for heavy target nuclei (p+^^^Ta 
at 4 GeV/c [3]), but even for light ones (p+^^Ne at 608 MeV/c |2J). The direct mecha- 
nism pp — i- AA can only explain less than 20% of the measured A production cross section 
for p+^^^Ta at 4 GeV/c [J]. Thus, the hyperon production in p annihilation on nuclei is 
dominated by the interactions of secondary particles produced in pN annihilation with the 
nuclear medium. Extensive theoretical calculations within the intranuclear cascade (INC) 
model Sj tend to overestimate the yields of both particle species, A and K^, and underesti- 
mate the yield ratio A/if^. The latter might indicate the in-medium enhanced cross sections 
of the strangeness exchange reactions, KN — )■ Yti {Y stands for a A- or S-hyperon). The 
enhanced hyperon production has been also interpreted in terms of p annihilations on the 
clusters of nucleons %m\ or even the formation of a cold quark-gluon plasma (QGP) [7|. 

Even more interesting is the double-strange (S=-2) hyperon production in antiproton 
annihilation on nuclei which has not yet been studied experimentally. The calculations 

[9|, [lOj, take into account only the direct mechanism A^A^ — > ^^ which can only account 
for a few percent of an inclusive S yield. However, the inclusive production of the S=-2 
hyperons in p annihilation on nuclei is very interesting by itself. This process should be a 
quite sensitive test for the unusual mechanisms of the p annihilation on nuclei, like the QGP 
formation, since it requires the simultaneous production of two ss pairs. 

In this work, we study strangeness production in p-nucleus collisions at piab = 0.2 — 15 

GeV/c and in p-nucleus annihilation at rest within the microscopic GiBUU tr ansp ort model 

ll| . We compare our calculations with all available experimental data [iMsl, ll2| on A, K'g 



and A production for in-fiight annihilation of p on nuclei and with earlier calculations within 
the INC models p, ll3|, [ij]. We also analyse the selected data set of ref. la] on vr^, proton 
and A production in p annihilation at rest on ^^N. In contrast to the INC models, the GiBUU 
model includes self consistent relativistic mean fields acting on baryons, antibaryons, kaons 
and antikaons. The selfconsistency means that the mean fields depend on the actual particle 



densities and currents. This has never been done in the previous transport calculations of 
antiproton- nucleus reactions. The selfconsistent potential fields are very important, for 
example, for a realistic treatment of annihilation on light nuclei, when a nucleus get almost 
destroyed by pions produced in annihilation, and the outgoing particles propagate in much 
weaker potential fields. 

We study in detail the mechanisms of strangeness production by decomposing spectra and 
reaction rates into the partial contributions from various elementary processes. Estimates for 
the production probabilities of hypernuclei are also given. Finally, we present the predictions 
of our model for the S hyperon spectra. We argue that the latter can be used to disentangle 
the hadronic and hypothetic QGP mechanisms of p annihilation on nuclei. 

The structure of the paper is as follows. Sec. Hi] contains a brief description of the GiBUU 
model with an accent on its new and/or improved ingredients, such as, e.g., strangeness pro- 
duction channels in NN collisions. In sec. Illlt we present the results of our calculations 
for the A, Kg and A production and compare them with experimental data and INC cal- 
culations. Then, we give our predictions for the S hyperon production at various p beam 
momenta. Finally, sec. II VI summarizes our work. 

II. MODEL 

The GiBUU model [ill] is a transport-theoretical framework which allows to describe a 
wide range of photon-, lepton-, hadron-, and nucleus-nucleus reactions. Below we concen- 
trate mostly on the model ingredients which govern strangeness production in antiproton- 
induced reactions. For other model details relevant for the present study we refer the reader 



to refs. 
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16|-ll9|. 



The GiBUU model solves the coupled system of kinetic equations for the different 
hadronic species i = N,A, Y, Y*, S, vr, K etc and respective antiparticles: 



1 
Po 






dxf^ \ !^ « I Q^ I Qp 



*fl 



/•,(x,p*) =/,[{/}], (1) 



where /j(x, p*) is the phase-space distribution function, p* = p — Vi is the kinetic four- 
momentum, Ff^'^ = d'^V^ — d'^Vf is the field tensor, and m* = irti + Si is the effective mass. 
The scalar potential Si = Qaicr is expressed in terms of the isoscalar a meson {J^ = 0~^) field. 
The vector potential Vi includes the contributions from the isoscalar u meson (J^ = 1~) 



field, isovector p meson (J^ = 1 ) field and the electromagnetic field A: 

Vi = g^iU + gpiT^p^ + qA . (2) 

Here, the isovector term oc p^ is included only for nucleons and antinucleons. The collision 
term /i[{/}] i^ ^^^ r.h.s. of Eq. ([1]) represents the contribution of binary collisions and 
resonance decays to the partial time derivative -^fi{x,'p*)- The term E^V/j(x, p*) in the 
l.h.s. of Eq. (Il]) is a usual drift term. The term oc -Q%ifi{x,p*) describes the deviation 
of particle trajectories from straight lines as well as acceleration/deceleration due to the 
meson mean fields. Without this term, the GiBUU model is basically reduced to a usual 
cascade model. The meson mean fields a, u, p and the electromagnetic field A are calculated 
from the corresponding Lagrange equations of motion neglecting time derivatives 16|, Il9 |. 
This calculation is selconsistent in the sense that the fields are induced by actual particle 
densities and currents which serve as source terms in the Lagrange equations. The full 
system of transport equations (II]) and meson field equations admits the energy-momentum 
conservation (c.f. ref. |l6|). 

In order to determine the meson-nucleon coupling constants and the self-interaction pa- 
rameters of the a field, we use the relativistic mean field (RMF) model in the NL3 version 
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20l |. The coupling constants of mesons with other nonstrange baryonic resonances (A, A^* 
etc) are set equal to the respective meson-nucleon coupling constants. For the meson-hyperon 
and meson-kaon couphng constants, we apply a simple light quark counting rule by putting 

(3) 
The coupling constants with the corresponding antiparticles are obtained as follows: 

g^B _ gcrY g^N /^n 

~;7^ — 'ly — gcrS — gcrK — ~~y ) l<3j 

guB _ g^jY _ g^jN ,„■. 

-177- — —7- — gujE — gujK — :r- 5 v"J 

gpN = igpN ■ (7) 

Here, "-B" denotes any nonstrange baryon, i.e. B = N,A, N* etc. The relative signs in 
Eqs. ©-([T]) are obtained from the G-parities of the meson fields. The factor C, = 0.22 is 
introduced in order to obtain the Schrodinger equivalent potential f/^ = —150 MeV for 



an antiproton at the zero kinetic energy and normal nuclear density po = 0.148 fm ^, in 



agreement with p- nucleus scattering phenomenology (c.f. 18| and refs. therein). One should 
keep in mind, however, that in experiment only the nuclear surface is tested due to a large 
annihilation cross section. Thus, an empirical information at po can only be obtained by 
model dependent extrapolations. 

In Table HI we collect the Schrodinger equivalent potentials of different particles evaluated 
by using the relations ©-(IT]) and the nucleon scalar {Sn = —380 MeV) and vector (V^ = 308 



MeV) potentials in nuclear matter at po 18|. The potential depths for nucleon, A hyperon 



TABLE I: The Schrodinger equivalent potentials of different particles at zero kinetic energy, 
Ui = Si + V^ + {Sf - {Vl^f)/2mi (in MeV), in nuclear matter at po- 
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-150 


-449 


-449 


-227 


-18 


-224 



and antikaon are consistent with phenomenology 2l|. The S potential is attractive, which 



contradicts the analysis of S^ atoms 21|. A weak attraction for kaons is also not supported 
by the analysis of the kaon flow from heavy ion collisions (c.f. |22[ and refs. therein), where a 
weak repulsion has been found. These drawbacks are the consequences of a simple treatment 
of particle potentials based on the same RMF model. We do not expect that they sensitively 
influence our results since the multiplicity of S hyperons is considerably smaller than the 
multiplicity of A hyperons while the kaon potential is weak anyway. The potentials of the 
S hyperons and antihyperons are still not restricted by any experimental data. 

The nucleus is modeled by employing a local density approximation. The momenta of 
nucleons are sampled uniformly within the spheres of radii 

PF,i = (37rVj)^^^ , i = p,n , (8) 



by using a Monte Carlo method. The density profiles of protons and neutrons 



£i 



V) are 



obtained from a selfconsistent solution of the relativistic Thomas- Fermi equations [19|. This 
makes the nucleus stable on the time scale of the order of several 100 fm/c, enough for the 
most reactions with nuclear targets. The Fermi motion of nucleons results in the smearing 
of particle production thresholds, which is especially important for the reactions with low- 
energy projectiles. 



For the simulation of an p-nucleus collision, at t = an antiproton is placed at the 
distance 5 fm + nuclear radius from the nuclear center along the collision axis. Such a 
distant initialization of the antiproton is needed in order to take into account the change of 
its momentum and trajectory under the action of attractive nuclear and Coulomb potentials 



18|. 



For the simulation of annihilation at rest, the initial radial position of the antiproton is 



chosen according to the probability distribution (c.f. |23|, |24| ) 



dP = C\Rniir)\^pirydr . (9) 

where Rniir) is the radial wave function of the antiproton in the Coulomb atomic state 
with quantum numbers n and / < n — 1, p{r) is the nucleus density profile, and C is a 



normalization constant. Quoting ref. 2J], a cascade of electromagnetic deexcitation of the 
p-atom, essentially through the {n,l = n — 1) states, emits X-rays, which permits us to trace 
the p down to the level where annihilation takes place. 

Once a nucleus and an antiproton are initialized, the system of kinetic equations ([1]) 
supplemented by the meson field equations is solved by the test particle method in the 



parallel ensemble mode 25|. Between the two-body collisions, the test particle centroids 



(r, p*) propagate according to the Hamiltonian-like equations (c.f. |18l. Il9j). The two-body 
collisions are treated in a geometrical picture: the test particles a and b from the same 
parallel ensemble are allowed to collide if they pass their minimal distance d^i^ during a 



given time step. Here (imin < yo'^^t/'^, where a^^^. is the total interaction cross section 
depending on the types of colliding particles and invariant energy ^/s. In the present work, 
we use vacuum cross sections neglecting their possible in-medium modifications. The final 
state of a two-body collision ab — )■ cde... is chosen by a Monte Carlo method according to 
the partial cross sections of various outgoing channels. The final state probability includes 
the products of Pauli blocking factors (1 — /) for the outgoing nucleons. These factors are 
calculated from the actual time-dependent neutron and proton phase-space distributions. 
Collisions between the secondary and primary particles as well as between two secondaries 
are taken into account. The mean fields are recomputed on every time step according to the 
modified particle densities and currents. 

In the case of an antinucleon-nucleon collision, the annihilation may result in a large 
number of various mesonic final states. This makes the direct parameterization of all partial 



cross sections practically impossible. Thus, we rely on the statistical annihilation model with 



SU(3) symmetry [26|,l27|. This model has been successfully applied in the calculations of pion 



and proton spectra from p annihilation on nuclei at 608 MeV/c within the GiBUU framework 



18| . According to the model 26|, l27|], the probability of the NN annihilation to a given 
final state meson configuration, which may include up to n = 6 mesons vr, rj, u, p, K, K, K* 
and K*, is defined as 

Wniy/s,Ii,...,In,Yi,...,Yn) = U;„(v^, Ji, ...,/„, Fi, ..., Fn) 

X a-a:;''a-a;''a7+"^a^-*+"^'* , (10) 

where /i,...,/„ and Yi,...,Yn are, respectively, the isospins and hypercharges of outgoing 
mesons, and Ui are the multiplicities of mesons of each type (z = •n^rj^oj etc). The quantity 
Wn is proportional to the phase space volume of a given final state and is calculated assuming 
that the incoming and outgoing hadrons can be exactly classified according to the SU(3) 
symmetry. The dimensionless parameters a^r, a^, ..., a/^-* break the exact SU(3) symmetry. 
They approximate the unknown parts of matrix elements and depend on the types of the 
particles and on their internal structure. For the annihilation channels without strangeness, 



the values a^r = 1,0,, = 0.13, a^^ = 0.18 and Op = 0.24 were determined in |28|] from the best 
agreement with the data on pp annihilation at rest and in flight at piab < 10 GeV/c J29|. 
For the parameters related to the strange mesons, we apply the beam momentum dependent 
expressions obtained from the fit of the pp — )■ K^X and pn — )■ K^K^Tf^ir" cross sections 
(see Fig. [T]for the pp case): 

a/r{/^.) = 0.07 (0.05) C(piab) , (11) 

where 

C4.2exp(-1.4piab) + 0.5exp(-0.16piab) + 1, ioi pp 
C(piab) = { , . , . _ (12) 

I 2.0exp(-0.7piab) + 0.5exp(-0.16piab) + 1, for pn 

with piab being the beam momentum (in GeV/c). The statistical annihilation model works, 
strictly speaking, only at high beam momenta, when the particle multiplicities are large. At 
low beam momenta, this model has to be supplemented by the phenomenological branching 
ratios of the different annihilation channels. For the channels without strange particles, this 

In the present work, we have extended the tables of 



has been already done in |26 

probabilities for various pp and pn annihilation channels at rest [26|, |27[ by including the 



channels with strange particles KK^ K*K + c.c, and K*K* (see Appendix |A]). In order to 
have a smooth transition from these empirical branching ratios to the description according 
to the statistical model as the beam energy grows, we determine by a Monte-Carlo method 
whether the statistical model itself or the empirical branching ratios are used to simulate a 
given NN annihilation event. The probability to choose the tables is 

Pat rest = max f 0, 1 - /lr^7 ] , (13) 

where ^Smax = 2.6 GeV is the maximum invariant energy up to which the annihilation tables 
at rest still can be selected (respective beam momentum piab = 2.5 GeV/c). At the invariant 
energies above ^Smax, the statistical model is used directly. The momenta of outgoing 
mesons in an NN annihilation event are always sampled microcanonically according to 
the available phase-space volume, regardless of whether annihilation tables at rest or the 
statistical model are applied to choose the flavours and charges of the outgoing mesons. 
This ensures a smooth behavior of the kinematics of produced mesons with increasing beam 
energy. 

Apart from the NN annihilation to mesons, we also include elastic (+charge exchange) 
scattering NN — )■ NN , A resonance production NN — t- iVA (+c.c.), and hyperon produc- 
tion NN — 7- AA, SA(+c.c), SS channels. The corresponding total and angular differential 
cross section parameterizations are obtained from the fits of empirical data. The cross 
sections are described in detail in Appendix [Bl 

At y/s > 2.37 GeV, the inelastic production in NN collisions is simulated with a 



help of the Fritiof model [30|. Exception are the processes NN — )■ NN, NN — )■ 
AA, SA(+c.c), SS which are either not included or not described well in Fritiof. Thus, 
we treat these processes separately according to their partial cross sections at any ^/s. 

Figure [H shows the inclusive cross sections of Kg and hyperon production in pp collisions. 
The partial cross sections of Kg production in the annihilation and nonannihilation events 
are also shown in the same figure. The experimental data on Kg production are described 
very well, including the nonannihilation channels, which are simulated by the Fritiof 
model. The hyperon production cross section is somewhat underestimated at the beam 
momenta ~ 4 GeV/c. This is mainly due to still underestimated cross section of the pp — >■ 
AAX inclusive channel above 3 GeV/c, which is practically missed in the Fritiof model. 
However, we find that, overall, the Fritiof model describes the inelastic production in pp 
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FIG. 1: Inclusive cross section of Kg production (upper left panel), the cross section of Kg produc- 
tion in annihilation and nonannihilation channels (lower left and upper right panels, respectively), 
and the inclusive cross section of the hyperon production (lower right panel) in antiproton-proton 



collisions as a function of the beam momentum. The experimental data are taken from refs. 



31 



2M- 



collisions better, than, e.g., the Pythia model [33|. Generally, the latter is successfully 
employed in GiBUU for the description of baryon-baryon and meson-baryon collisions at 
high invariant energies {^/s > 2.6 GeV and ^/s > 2.2 GeV, respectively). We have to only 
admit one problem. The Pythia model does not include K^ and K^ in its list of the 
possible incoming particles. Thus, both of them are replaced by K^ in GiBUU every time 
when Pythia is used. This leads to ~ 5% violation of the total strangeness conservation in 
our calculations, which, however, is accurate enough for the present exploratory studies. 



The processes of a hyperon scattering, AA^ -)■ AA^, AA^ ^ 'EN, 'EN -^ 'EN, EN -)■ EN, 
EN —7- AA, EN — 7- AS, and of the strangeness exchange on nucleons, KN — > Yn, KN — )■ 
EK, are taken into account in the modeL When present, their empirical cross sections 
have been suit ably parameterized or the fits to the existing theoretical calculations have 



been done 



34. 



35( 1 ■ As we will see later on, the strangeness exchange processes are very 
important for the hyperon production in antiproton-induced reactions on nuclei. In GiBUU, 
the strangeness exchange reactions KN — ?► Ytt are partly mediated by the 5* = — 1 hyperon 



resonance formations and decays KN — > F* — ;■ Yn [35 1. 

The associated hyperon production is included in GiBUU via reaction channels tiN — >■ 
YK, r]N — )■ YK, pN — )■ YK and uN — )■ YK. The cross sections vrA^ — )■ YK have been 



parameterized according to Ref. |36j. The ry-induced associated hyperon production cross 
sections on proton have been reconstructed from the respective cross sections of the 7r°- 
induced processes at the same invariant energy y/s by utilizing the detailed balance relations 



371] and isospin invariance: 



a{r]p -^ K+A) = a{7i'>p ^ K+A)^^ , (14) 

Pr}N 

aivp ^ K+E^) = a{n^p -^ K+E°)^^ , (15) 

PriN 

(t{7]p -^ K°S+) = 2a{7]p -> fs:+S°), (16) 

where Pj^n and Pr^N are the center-of-mass (cm.) momenta for the corresponding initial 
channels calculated at the same ^/s. Similar formulas have also been used for the up initial 
channel. For the p-induced reactions on proton, we have assumed 

a{pp -^ KA) = ainp -^ KA)^^ , (17) 

PpN 

a{pp -^ KE) = a{np -^ KE)^^ , (18) 

PpN 

where the isospin states of all particles match each other. The cross sections of the t]-, uj- 
and p-induced reactions on neutron have been obtained by using the isospin reflection from 
the corresponding cross sections on proton. 

We note, flnally, that the bubble chamber data on A production contain also the admix- 
ture of A's produced by the decays S" — )■ A7. These decays are not included in GiBUU, 
since the S° life time r = 7.4 x 10~^° s ~ 222000 fm/c is much longer than the typical 
hadron-nucleus reaction time scale of ~ 100 fm/c. Thus, in the present calculations, we 
simply add the S° yield to the A yield. 
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III. RESULTS 

A. Time evolution of hyperon production 

In our calculations, baryons experience the action of attractive mean field potentials. 
Slow hyperons get captured inside the residual excited nuclear system. This system may 
evaporate particles and/or decay into fragments, some of which will be single- or double-A 
hypernuclei. It is natural to assume that the fragmentation and evaporation will not change 
much the total yield of hypernuclei; they may, however, affect the production of a given 
hypernuclear species. 

In order to distinguish the hyperons outside and inside the residual nucleus, we have 



applied a simple criterion based on the relative distance between particles [38|: a particle 
is outside the nucleus if it is separated by the distance larger than some critical distance dc 
from all other particles of the nucleus. Otherwise, the particle is inside the nucleus. Provided 
the evolution time is long enough, the result should not be much influenced by the choice 
of dci if the latter is larger than the internucleon spacing ~ 1 — 2 fm. This is illustrated 
in Fig. [2] where we show the number of (A + S°) hyperons outside and inside the residual 
nucleus per one annihilation event as a function of time for p(225 MeV/c)+^^C collisions 
for various choices of the critical distance. After ~ 200 fm/c the both numbers change very 
slowly indicating the presence of really captured A hyperons inside the attractive potential 
well. For comparison, we also present the results for the bound and unbound hyperons in 
the same figure. In this case a hyperon Y is considered to be bound (unbound) if Ey < my 



{Ey > my), where Ey = a/(p*)^ + {rriy)^ -\- Vy is its single-particle energy. This criterion 
allows to identify the captured hyperons somewhat earlier. However, after 200 fm/c it is 
very close to the criterion according to d^ = 3 fm. 

Of course, the capture may only happen if a particle experiences the action of an attractive 
potential. Fig. |3] compares the time dependence of the number of A's inside and outside the 
nucleus for calculations with and without A-potential. As expected, in the calculations 
without A-potential, the number of hyperons inside the nucleus quickly drops with time 
indicating that there are no captured A's in this case. 

In the following, we always fix rfc = 3 fm as in ref. 



and calculate the time evolution 



until 200 fm/c. Our results for the number of free hyperons, i.e. those outside the residual 
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FIG. 2: The number of (A + S°) hyperons outside (out) and inside (in) the residual nucleus per 
annihilation event as a function of time for p annihilations at 225 MeV/c on ^^C. The calculations 
are done for three different critical distances dc and for the criterion according to the hyperon 
binding energy. The results are weighted with the impact parameter of an incoming p. 

nucleus, change only by ~ 10% if we further increase the evolution time (c.f. Fig. |2]). As 
demonstrated in Fig. HJ this change concerns only slow hyperons, while the yields of fast 
particles are practically stable. In particular, the kaon yields and spectra are not influenced 
by any further increase of the evolution time. 



B. Annihilation at rest 

We selected the ASTERIX@LEAR data [151] on charged pion, proton and A production 
from p annihilation at rest on ^^N. According to ref . [39| , the last observable transition in 
light antiprotonic atoms is 4 — )■ 3. Thus we assumed that the antiproton occupies mainly 
the n = 3,1 = 2 level immediately before annihilation which we used as an input in our 
calculations (see Eq. ([9])). We have checked, however, that our results are changed by 
only few percent if the quantum numbers n = 4,1 = 3 are chosen for the antiproton wave 
function. 
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FIG. 3: Same as in Fig. [2] but for calculations with and without A-potential at the fixed critical 
distance dc = 3 fm. 

Fig. [5] shows the vr^, p and A kinetic energy spectra in comparison with our calcula- 
tions. In ref. 15|], the measured spectra were fitted by the relativistic Maxwell-Boltzmann 
distribution 

(19) 



E^r^ = Aexp{~E/Eo) 



p^dp 

with the normalization A and the temperature Eq being fit parameters. We observe that the 
high energy parts of calculated spectra — except for very high energy (-Ekin > 0.6 GeV) pion 
spectrum — agree with the data and with the Maxwell-Boltzmann formula ( TTOjl reasonably 
well. A slight underprediction of the high-energy pion spectrum can be traced back to the 
case of pp annihilation annihilation at rest (c.f. Fig. B.52 in ref. [llj). More significant 
is the deviation from the data at small kinetic energies. We predict the strongly enhanced 
evaporative emission of the low-energetic particles, in-particular, p and A, in contrast with 
the ASTERIX data. The enhanced emission of low-energy nucleons is also present in the 
INC calculations of ref. 2J] for the p annihilations at rest on ^^Mo. 

The calculated pion spectrum also has a clear two-component structure, which seems to be 
absent in the pion spectra measured by ASTERIX. The high energy pions (-Ekin > 300 MeV 
or piab > 400 MeV/c) originate directly from NN annihilation almost without rescattering 
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FIG. 4: Kinetic energy spectra of (A + S^) hyperons (left panels) and K'g mesons (right panels). 
The spectra of free particles selected according to the two different criteria, Ei > rrii (i = Y, K, 
upper panels) and dc = 3 fm (middle panels), are plotted at different times. Also the comparison 
of the spectra obtained with these two criteria with the spectrum of all particles, i.e., without any 
restrictions is shown at 200 fm/c (lower panels). The number of Kg mesons has been determined 
from the numbers of K'^ and K'^ as A^^o = ^{Nj^o + Nj^o). 
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FIG. 5: Kinetic energy spectra of negative pions, protons and (A + S^) hyperons for the p annihila- 
tion at rest on ^^N. The GiBUU calculations are performed for the n=3, 1=2 and n=4, 1=3 states of 
the antiproton. Experimental data and their fits according to the relativistic Maxwell-Boltzmann 

distribution of Eq. (J19p are from ref. [l5|. The slope temperatures Eq for the experimental spectra 
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are also given on the plots. 
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on nucleons. The low energy pions are mostly the products of the vrA^ — )■ A — ;■ ttN processes. 



This structure is present in the CALLIOPEQlLEAR data [40|] for the pion momentum spectra 
from 608 MeV/c antiproton annihilation on ^^C and ^ss^j ^j-^jch agree with the GiBUU 



calculations very well 



18|. 



Although the radial distributions of annihilation points are somewhat different for annihi- 
lation at rest and low-energy annihilation in-flight, the an gle- integrated momentum spectra 
of emitted particles are expected to change only a little [24]. Therefore, we attribute the 
above discrepancies to the reduced acceptance of the ASTERIX spectrometer for the low- 
momentum particles, as also mentioned by the authors themselves in ref. |l5l |. 
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C. Annihilation in-flight 

We start from the lowest beam momenta and consider the reactions p(0-450 
MeV/c)+^2Q^48rpj^i8irj^ g^^^ 208pb measured at BNL [1]. Figure E shows the calculated 

target mass dependence of the free A production probability per annihilation event together 
with the data. Note, that the ^-annihilation cross section (Tann ex A^'^, i.e., it grows with 
the target mass number. However, we got rid of this enhancement by dividing out (Tann from 
the A production cross sections. We also present the calculated production probabilities of 
nuclear systems with one and two captured A hyperons in the same reactions. In calcu- 
lations, the beam momentum was fixed at 225 MeV/c. The free A production probability 
is weakly sensitive to the target mass number and agrees with experiment. However, the 
probabilities of one and two captured A production grow with the target mass number by 
almost one order of magnitude from the lightest (^^C) to the heaviest (^°*Pb) target. This 
is expected since, in heavier targets, the produced hyperons are more efficiently decelerated 
by rescattering on nucleons. The detailed calculations of hyperfragment production in p 
annihilation on nuclei within the coupled GiBUU -t- statistical multifragmentation models 



are presented in Ref. |34|. 



Figure [7] shows the rapidity distributions of A hyperons and Kg mesons from p+^'^Ne 
collisions at 608 MeV/c in comparison to the data and to the INC calculations from |5|. 
As we observe, both models agree within ~ 30%. The experimental A rapidity distribution 
is rather well described, while the calculated Kg spectra are clearly above the data. Since 
the direct channel of hyperon production, pp — > AA {p\^l = 1.439 GeV/c), is closed at 
Piab = 608 MeV/c, the hyperons can only be produced in strangeness exchange processes 
K{K*)N — 7- Y{Y*)tt, in hyperon resonance formation reactions KN — )■ Y*, or in meson- 
nucleon collisions {uN — )■ YK mostly). Our calculation produces somewhat more Kg's than 
the INC calculation [5|] does. This is partly due to taking into account the target destruction 
by p annihilation in our calculations, which reduces the chances for an K to be absorbed. 
There is also another possible reason for the different results. Both models contain the ^A^ 
elementary cross sections fitted to the experimental data. However, the important difference 
is that in GiBUU the Y* resonances have finite life times ~ l/F and, therefore, can escape 
out of the nucleus. This effectively reduces the K absorption since a Y* excited as the 
intermediate state of KN scattering, KN -^ Y* ^ KN, "hides" an K from interactions 
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FIG. 7: The summed rapidity distributions of A and S*^ hyperons (upper panel) and the rapidity 
distribution of K^ mesons (lower panel) produced in p annihilations on ^^Ne at 607 MeV/c — 
thin solid lines. The integrated distributions give the number of particles per annihilation event. 
The partial contributions to the hyperon rapidity spectrum ("M" means any nonstrange meson 
vr, 77, p, uj) are depicted by lines with open symbols. The INC calculations from ref. a] are shown 
by dashed lines. The experimental data (histograms) are taken from ref. [4]. The cm. rapidity of 
an antiproton and a nucleon ypjy = 0.305 is shown by arrows. 
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with other nucleons. 

The (A + S'') rapidity spectrum shown in Fig. [7] is decomposed to its partial components 
according to various elementary production channels of A or S°. We see that the largest 
contribution is given by the hyperon rescattering on nucleons with flavour and/or charge 
exchange, in-particular, by the exothermic processes S+n — > Kp and S^p — )■ An. The 
main channel of the hyperon production, however, is the strangeness exchange processes 
K{K*)N -> F(r*)7r. The rescattering A(So)A^ -^ A(So)A^ somewhat distorts the shapes 
of these partial contributions by shifting the produced A's(S°'s) towards smaller absolute 
values of rapidity (see also Fig. [3 below). 

A more detailed information on the relative importance of various hyperon production 
channels is given in Table [TTl As one can see by inspecting Table [TTl ~ 80% of the total 
Y and Y* production rate (without Y and Y* absorption contributions) is caused by the 
KB —7- YX, KB — > Y* and KB — t- Y*tt secondary processes. Nevertheless, the associated 
production MB — )■ YK constitutes the remaining 20% and, therefore, is quite important 
as well. In contrast, 80% of the kaon and antikaon production rate is caused by the direct 
mechanism BB — )■ KKX. The INC calculation gives larger contributions of the antikaon 
absorption KB — t- YX and associated production MB — t- YK. At the same time, the 
hyperonic decays of Y* resonances present in GiBUU counterbalance the smaller A and S 
hyperon production by other channels. This leads to rather similar final results for the both 
models. 

It is quite interesting to observe from Table Hi] that the irB — ;■ YK channels make rela- 
tively small contributions to the A and S production rates with respect to the uB — )■ YK 
channels even though there is an abundant pion production in p annihilation on nuclei. The 
reason is that in p annihilation at rest or at low p beam momentum most of pions are pro- 



duced with momenta ~ 0.2 — 0.4 GeV/c (c.f. Fig. 2 in ref. 18|, where the 7r+ momentum 
spectrum is shown for p+^'^C collisions at 608 MeV/c). This is well below the threshold 
pion momentum of 0.895 GeV/c for the ttN — > AK reaction on a nucleon at rest. On the 
other hand, an u meson is produced in a rather large fraction (~ 20%) of NN annihilation 
events at rest (c.f. ref. [27|). This favors the exothermic uN — )■ YK reactions, which have a 
large cross section at low u momentum. The situation is, however, different for the energetic 
p-nucleus collisions (see Table UTTl below) . 

At 4 GeV/c, the agreement of our calculations with the experimental Kg spectrum be- 

19 



10' 



^ 10' 

T3 
T3 

10^ 



10' 



p(4 GeV/c)+^^^Ta->(A+I°)X 



#•• 



\* i^pN 




2.0 



p 

T3 



10' 



10' 



10' 



10" 




-2.0 



-1.0 



0.0 



1.0 



2.0 



FIG. 8: The summed rapidity distributions of A and T,^ (upper panel) and the rapidity distribution 
of Kg (lower panel) for the collisions p(4 GeV/c)+"'^^"'^Ta are shown by solid lines. The partial 
contributions to the hyperon rapidity spectrum ("M" means any nonstrange meson 7r,r],p,uj, "B" 
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ref. f3[. The antiproton-nucleon cm. rapidity yp^ = 1.078 is shown by vertical arrows. 
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TABLE II: Partial contributions of the different reactions to the multipHcities of A, T, K and K 
per 1000 p absorption events on ■^'^Ne at 608 MeV/c. In the calculation of reaction rates we did 
not distinguish kaons from K*^s and antikaons from ^*'s. The negative numbers correspond to 
the absorption reaction channels for a given particle. Results of the INC calculations [5] are given 
in brackets. 



BB ^ AK (0) BB^llK (0) BB -^ KX 57.7 (50.0) BB -^ KX 57.7 (50.0) 

ttB -^ AK 0.6(1.4) ttB-^T.K 0.6(1.0) T.M -^ KX 0.2 (p-^KK 1.9 

r]B -^ AK 0.2 (1.5) r]B -^ T.K 0.3 (0.4) Y*M -^ KX 0.5 ttB -^ AK 0.6 (1.4) 

pB -^ AK 0.9 pB-^T.K 1.5 4>-^KK 1.9 rjB -^ AK 0.2(1.5) 

ojB -^ AK 1.5 (5.0) ojB -^ SK 3.1 (4.0) MM -^ KK 3.2 pB -^ AK 0.9 

KB -^ AX 8.4 (10.0) KB -^ T.X 10.7 (16.0) KK -^ MM -1.2 uB -^ AK 1.5 (5.0) 

Y.B -^ AX 9.9 (10.0) Y.B -^ AX -9.9 (-10.0) KB -^ AX -8.4 (-10.0) ttB -^ SK 0.6 (1.0) 

Y*B -^ AX 1.3 Y*B -^ T.X 1.4 KB -^ T.X -10.7 (-16.0) r]B -^ T.K 0.3 (0.4) 

AB -^ T.X -0.9 (-0.9) AB -^ T.X 0.9 (0.9) KB -^ Y* -17.0 pB -^ T.K 1.5 

AM -^ Y*X -1.5 SM -^ Y*X -1.5 Y* -^ KN 6.4 uB -^ T,K 3.1 (4.0) 

Y* -^ Anir]) 5.9 Y* -^ Svr 7.6 KB -^ Y*-k -2.7 MM -^ KK 3.2 

BB -^ AX 0.2 EM -^ KX -0.2 KK -^ MM -1.2 

BB -^ KX 0.3 

Sum A 26.5 (27.0) Sum S 14.5 (12.3) Sum K 29.9 (24.0) Sum K 70.6 (63.3) 

comes better, as one can see from Fig. [HI However, the (A -|- S°) yield is now underpredicted 
at intermediate rapidities y ~ O.bypiy, where ypjy is the rapidity of an antiproton-nucleon 
cm. system. From the spectrum decomposition into partial channels, which is also given in 
Fig. IHl we observe qualitatively the same picture as at lower beam momentum 607 MeV/c 
(Fig. [7] above). Some changes with respect to the lower beam momentum are, however, visi- 
ble: The hyperon resonance Y* decays grow to the second place in importance. The relative 
contribution of nonstrange meson-baryon collisions is increased and becomes as important 
as that of strangeness-exchange reactions. Also a nonnegligible contribution of direct BB 
channels appears now. 

The disagreement with the experimental rapidity spectrum of (A-|-E°) hyperons at y ^ 0.5 
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might be due to still underpredicted contribution of the Y* — )■ Yn decays. Indeed, this 
contribution has a broad maximum at y ^ 0.3, i.e., close to the maximum of the measured 
rapidity spectrum. Since our calculations tend to overestimate the Kg production, we can 
assume that the KN — )■ Y* cross sections should be larger. In-particular, the total cross 
section K^n — > X (or the same cross section of the isosp in-reflected channel K^p — )■ X) is 



underestimated at y/s < 1.7 GeV (see Fig. 2.15 in ref. |35|). 

Another reason for the deviation with the experimental rapidity spectrum of (A + E*^) 
hyperons is related to rather uncertain YN — )■ YN cross sections. Especially the AA^ — )■ AA^ 

n 

cross section is important, as has been noticed, first, by Gibbs and Kruk in ref. |14j |. 



These authors have applied their INC code 13|, [l^ to the reaction p(4 GeV/c)+^'^^Ta and 
reproduced the measured (A + S°) and Kg rapidity spectra very well. They assigned a 
constant elastic hyperon-nucleon cross section of 14 mb, which is about two times less than 
the Ap elastic cross section at piab — 0.5 GeV/c (corresponding to the peak position y ~ 0.4 
of the measured A rapidity spectrum, c.f. Fig. [8]) in the parameterization of Cugnon et al. [5|. 
Since we have also applied the latter parameterization in the present calculations, this largely 
explains the shift of our (A + S*^) rapidity spectrum to smaller rapidities with respect to the 
measured spectrum. Fig. [9] demonstrates the sensitivity of our calculations to the hyperon- 
nucleon cross sections. The calculation without hyperon-nucleon rescattering produces the 
peak position of the calculated (A + S°) rapidity spectrum close to the experimental one. 
However, the spectrum at large forward rapidities is now overestimated. 

The stopping power of the nuclear medium with respect to the moving hyperon depends 
not only on the integrated elastic hyperon-nucleon cross section, but also on its angular 
dependence. For simplicity, we have chosen the hyperon-nucleon cross sections to be isotropic 
in the cm. frame. This also contributes to the somewhat too large deceleration of the 
hyperons. On the other hand, the total yield of (A + S*^) hyperons can be enhanced with 
increased charge exchange S^p — )■ An and S+n — t- Ap cross sections which are rather poorly 
known from experiment. 

Detailed information on the various production rates at 4 GeV/c is collected in Ta- 
ble mil One can see from Table IIIII that with a relative contribution of ~ 72% the kaon 
absorption reactions KB — t- YX, KB -^ Y* and KB — )■ Y*tt are still dominating in the 
total Y and Y* production rate, which is somewhat smaller relative contribution than in 
the case of p{607 MeV/c)^°Ne. This is expected, since at higher p beam momenta more 
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FIG. 9: Rapidity distribution of (A + S*^) for antiproton collisions at 4 GeV/c with ^^^Ta calcu- 
lated with (solid line) and without (dotted line) taking into account hyperon-nucleon rescattering. 
Experimental data are from [3|]. 

strangeness production channels open. The nonstrange meson-baryon colUsions provide the 
second largest contribution of ~ 23% to the Y and Y* production rate. The hyperon 
production in antibaryon-baryon collisions (including direct channel) and in baryon-baryon 
collisions contributes only ~ 3% and 2%, respectively, to the same rate. On the other hand, 
as before, in the case of p{607 MeV/c)^''Ne, the antibaryon-baryon collisions dominate in 
the total K{K) production rate contributing ~ 60%. It is interesting, that the meson-meson 
reactions MM -^ KK are rather important and contribute ~ 20% to the K and K pro- 
duction rates. This means that the KK pair production processes in mesonic cloud created 
after p annihilation should also be taken into account on equal footing with other secondary 
production channels. Note, however, that in GiBUU the two particles produced in the same 
two-body collision or resonance decay event are allowed to collide only after at least one 
of them collided with another particle not involved in this event. This is done in order to 
avoid multiple chain reactions between the correlated products of the same elementary event, 
which would otherwise lead to double counting in the production processes and, moreover, 
violate the molecular chaos assumption. 
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experimental data are taken from ref. [3]. 
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TABLE III: The same as Table [Ul but for reaction p{4 GeV/c)+^^^Ta. 



BB -^ AX 3.6 (13.0) BB -^ SX 2.2 (0) BB -^ KX 109.9 (178.0) BB -^ KX 110.1 (178.0) 

ttB -^ KK 6.6 (30.0) -kB -^ T,K 10.1 (36.0) SM -^ KX 1.1 </> -^ KK 2.8 

r]B -^ AK 1.2 (7.0) r]B -^ T,K 2.1 (4.0) Y*M -^ KX 3.8 ttB -^ AK 6.6 (30.0) 

pB -^ AK 2.5 pB -^ T.K 6.0 (p -^ KK 2.8 rjB -^ AK 1.2 (7.0) 

ujB -^ AK 2.0 (4.0) uB -^ T,K 7.5 (8.0) MM -^ KK 32.4 pB -^ AK 2.5 

KB -^ AX 20.1 (46.0) KB -^ SX 13.4 (60.0) KK -^ MM -4.5 coB -^ AK 2.0 (4.0) 

S5 ^ AX 47.5 (76.0) T,B -^ AX -47.5 (-76.0) KB -^ AX -20.1 (-46.0) ttB -^ T.K 10.1 (36.0) 

Y*B -^ AX 11.4 T,B -^ Y*X -1.3 KB -^ EX -13.4 (-60.0) r]B -^ T,K 2.1 (4.0) 

AB -^ SX -12.5 (-26.0) Y*B -^ EX 15.0 iffi -^ Y* -60.9 pS ^ J^K 6.0 

AS ^ y*X -2.5 AB -^ SX 12.5 (26.0) Y* -^ KN 22.4 coB -^ T,K 7.5 (8.0) 

AM ^ y*X -10.9 EM -^ Y*X -7.4 i?5 ^ Y*tt -25.0 MM ^ KK 32.4 

Y* -^ ATT{r]) 33.4 Y* -^ Stt 27.2 KK -^ MM -4.5 

55 ^ AX 2.4 SM ^ KX -1.1 YB -^ KX 3.8 (11.0) 

BB -^ SX 1.2 BB -^ KX 3.9 

Sum A 104.8 (150.0) Sum S 39.9 (58.0) Sum K 48.5 (72.0) Sum K 186.5 (278.0) 

Figure [TO] shows the inclusive transverse momentum spectra of (A-|-S°), Kg and (A-|-S°) 
from p(4 GeV/c)-|-^^^Ta collisions. The (A-(-S°) production at low pt is clearly enhanced due 
to large cross sections of the exothermic strangeness exchange processes KA^ — )■ Yn induced 
by slow antikaons. Moreover, the produced hyperons are decelerated by rescattering on 
nucleons. Our calculations reproduce the shapes of the experimental p^-spectra rather well. 
There is a remarkable agreement for the absolute values of a (A-|-S°) yield at high pf. Again, 
the yield of Kg is somewhat overestimated, while the yield of antihyperons is underestimated 
by our calculations. 

In Table HV] we summarize the results of comparison of our calculations with INC cal- 
culations and with experimental data at 0.6 and 4 GeV/c. Both models overestimate the 
strange quark production and underestimate the ratio o"a+so /<7_ft:o . 

To our knowledge, the latest measurement of neutral strange particle production from 
high-energy antiproton interactions with nuclei was performed at BNL using the Multipar- 
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TABLE IV: The cross sections of the (A + S*^) hyperon, K^ meson and strange quark production 
(mb) in comparison with the INC results from [5] and experimental data from [3, y] for reactions 
p(607 MeV/c)+^'^Ne and p(4 GeV/c)+^^^Ta. The "experimental" cross sections of the s-quark 
production are, actually, evaluated in [5(] based on some simple relations to the really measured 
cross sections. 



system cta+so o-ro o'A+so/o-a'O <^s 

p+20Ne GiBUU 18.7 19.2 1.0 50.0 

INC 19.6 13.7 1.4 39.7 

exp 12.3±2.8 5.4±1.1 2.3±0.7 16.9 

p+i8iTa GiBUU 154 121 1.3 387 

INC 275 142 1.9 454 

exp 193 ±13 82 ±6 2.4 ± 0.3 277 ±21 

tide Spectrometer (MPS) facility [12]. Fig. [11] shows the inclusive cross sections for A, 
Kg, A and strange quark production in collisions of antiprotons at 5, 7 and 9 GeV/c with 
carbon, copper and lead targets in comparison to the MPS data from [12|. Also the INC 
model [isl, llJ] results given in ref. pjj are shown in Fig. [HI The strange quark (or ss pair) 
production cross section has been calculated consistently with ref. [l2|, i.e. according to the 
approximate formula 

(^s = ^{4:(Tko +aA + CTso + o^A + c^Eo) ' (20) 

where a^o = {o-ro + <7i?o)/2. There is a fair overall agreement of GiBUU calculations with 
data. In particular, A and A production on the carbon target is described very well by 
GiBUU, while for the heavier targets we somewhat underpedict A and A production. The 
Kg production cross section on the carbon target is two times overpredicted by GiBUU. On 
heavier targets, the agreement with experiment on Kg production becomes better, but the 
slope of the beam momentum dependence of the if^-production cross section seems to be 
overpredicted. We note, that the data on the inclusive cross sections have been obtained 
by integration over rapidity region with good acceptance and extrapolated to the An solid 
angle using the INC calculations. This is partly responsible for a better agreement of the 
INC results with this experiment. Fig. [12] shows the rapidity distributions of A, Kg and 
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FIG. 11: Inclusive cross section of A, X^, A and strange quark (s) production for antiproton inter- 
actions at piab = 5.2, 7.0 and 8.8 GeV/c with ^^c, ^'^Cu and ^ospb targets. The INC calculations 
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FIG. 12: Rapidity distributions of A, Kg and A from p interactions with ^'^Cu at 8.8 GeV/c. The 
INC calculations and experimental data are from ref. [12| . 
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A from antiproton collisions with copper target at 9 GeV/c. As one observes, the GiBUU 
calculations agree with the data quite well, except the underprediction of A yield at y ~ 0.5. 
We also observe a rather close agreement of GiBUU and INC results, which means that the 
selected observables are, actually, not very sensitive to the model details, once the elementary 
cross sections are adjusted to the experimental input. 

D. 5 = — 2 hyperon production 

Figure US] shows the inclusive momentum spectrum of H hyperons from p+^^^Au collisions 
at 3 GeV/c together with the partial contributions from various S production channels. In 
performing this decomposition, we did not distinguish kaons from antikaons. As one can 
see from Fig. [131 the (anti)kaon-baryon collisions deliver the main contribution ~ 35% to 
the S production, mainly due to the double strangeness exchange channel KN — )■ KS. 
The decays H* — ;■ Hvr — especially important at low transverse momenta of S — make the 
second largest contribution ~ 26% to the S production. It is interesting, that the (anti)kaon- 
hyperon collisions, which are collisions between the secondary particles, contribute also quite 
appreciably, ~ 17%. Other reaction channels are of relatively minor importance for the 
inclusive S production. For example, the direct channel NN -^ SS contributes ~ 6% only; 
this channel is of primary importance for the planned PANDA experiment on the double-A 
hypernuclei production at FAIR |8|, |9|. 

Figure [H] shows the rapidity spectra of S~ hyperons together with the (A + H^) and Kg 
rapidity spectra from p+^^^Au collisions at 3, 9 and 15 GeV/c. The S~ spectra are about 
two orders of magnitude below those for the (A + S°) and Kg production. They are peaked 
at y ^ 0.5, 0.9 and 1.2 for the beam momenta of 3, 9 and 15 GeV/c, respectively. However, 
the (A + S°) spectra are always peaked near the target rapidity, y = 0, even at the largest 
beam momentum. This is because the hyperon production is dominated by the KN — )■ Yir 
processes with slow antikaons. Moreover, at 3 GeV/c, also the Kg spectrum has a broad 
maximum at the target rapidity. 

The experimental ia.ct that the (A + S°) and Kg rapidity spectra from p+^^^Ta collisions 
at 4 GeV/c are peaked at nearly the same rapidities (c.f. Fig. |8]) has been interpreted 
by Rafelski [7] as the manifestation of a common production source for strange particles, 
i.e. of an annihilation fireball with the baryon number ~ 10. The large baryon number is 
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FIG. 13: Inclusive momentum spectrum of H hyperons fromp(3 GeV/c)+^^^Au collisions (the solid 
line). The partial contributions according to the parent particles of a H hyperon are depicted by 
various lines with symbols. In key notations, "M" denotes any nonstrange meson TT,rj,p,uj, "B" — 
any nonstrange baryon N, A, N* etc.. The symbol '^K{K*)" denotes both kaon (K*) and antikaon 

{K*). 

due to absorbing nucleons by the propagating fireball. Rafelski assumed that the fireball 
undergoes the transition to the supercooled QGP state and then hadronizes. The rapidity 
spectra of S hyperons would also be peaked at the fireball rapidity if the fireball mechanism 
would dominate. In our model, a purely hadronic picture emerges instead, where the S 
production is dominated by the double strangeness exchange processes of KN — )■ KE type. 
The latter are endothermic and require the momentum of incoming antikaon in the rest 
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FIG. 14: Inclusive rapidity spectra of H^, (A + H,^) and K'^ from p+^^'^Au collisions at 3, 9 

and 15 GeV/c. The spectra of H^ are multiplied by a factor of 10. The pN cm. rapidity 

ypN = 0.940, 1.479 and 1.733 for pi^b = 3, 9 and 15 GeV/c, respectively, is shown by vertical 

arrows. 
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frame of a nucleon placed above the threshold value pthr — 1-05 GeV/c corresponding to 
the KN cm. rapidity rapidity of 0.55. This makes the S rapidity spectra shift to higher 
rapidities. In contrast, the S = —1 hyperon production is dominated by the exothermic 
strangeness exchange KN — )■ -kY . The cross section of this process grows with decreasing 
antikaon momentum in the nucleon rest frame. This favors the isotropic production of slow 
hyperons in the target nucleus rest frame. 

Finally, Fig. [15] presents the laboratory momentum spectra of the different strange par- 
ticles in the forward (0 = 0° — 60°) and transverse (G = 60° — 120°) directions for p+^^^Au 
collisions at 3, 9 and 15 GeV/c. The spectra of (A + TP) hyperons and K^ are rather close 
to each other both in forward and transverse directions. It is, moreover, interesting that the 
high-momentum slopes are similar for all considered particles. However, the production of 
the low-momentum S~ hyperons is suppressed due to the dominating S production via the 
double strangeness exchange in KN — )■ K'^X processes. 

IV. SUMMARY AND CONCLUSIONS 

To summarize, we have studied the strange particle production in p-induced reactions 



on nuclei by applying the GiBUU model. We have considered 



reactions and confronted our results with experimental data 



lim, 



3oth at-rest and in-fiight 



12 



15| on neutral strange 



particle yields and spectra. We have also compared our calculations with the earlier INC 
calculations of Cugnon-Deneye-Vandermeulen |5| and Gibbs-Kruk [l^. Finally, the model 
predictions for the S hyperon production are given. 

So far the main motivation for the experimental studies of strangeness production in 
antiproton-nucleus reactions has been to find the signatures of QGP production, which can 
be complementary to similar studies in high-energy heavy-ion collisions. Overall, our results 
are in a fair agreement with existing experimental data on K^, (A + TP) and (A + TP) 
production for p annihilation on nuclei and with INC models. There seems to be, indeed, 
no striking disagreements with data on inclusive A, K^ and A production at high beam 
momenta, which could point to the exotic mechanisms, as it has been already concluded in 
refs. 



12 



14- 



There are, however, some systematic deviations the origins of which remain to be better 
understood. The GiBUU and INC models are based on the same hadronic cascade picture 
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FIG. 15: Inclusive laboratory momentum spectra of H^, (A + S*^) and K^ from p(3, 9, 15 
GeV/c)+^^^Au collisions in the two ranges of the polar angle: = 0° — 60° and = 60° — 120°. 
The spectra of H^ are multiplied by a factor of 10. 
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but differ in several details like, e.g. the mean field potentials and the treatment of Y* 
resonances. There is a clear tendency to underestimate the A/ Kg ratio in our calculations. 
This tendency is present also in the results of INC calculations, although somewhat less 
pronounced. This problem might be caused by possible in-medium effects on the KN — >■ 
YX cross section. Poorly known hyperon-nucleon cross sections also strongly influence the 
observed hyperon spectra. The latter cross sections are also of more general interest, since 
they are related to the hyperon-nucleon interaction relevant in hypernuclear studies. To 
disentangle various possibilities, it would be very useful to measure the rapidity spectra of 
not only neutral strange particles, (A + S°) and K^, but also of the charged ones, S^ and 

The exotic mechanisms of strangeness production may, in-principle, manifest themselves 
more clearly in the 5* = —2 hyperon production. We have calculated the inclusive production 
of H hyperons from collisions p+^^^Au at 3-15 GeV/c. The H~ rapidity spectra are strongly 
shifted to forward rapidities in the laboratory system, since the S production via the domi- 
nating endothermic KN — )■ KS channel requires a high-momentum initial antikaon. This is 
in a contrast to the strangeness production mechanism from a moving annihilation fireball 
proposed by Rafelski [7]. Thus a simultaneous measurement of the single- and double-strange 
hyperon rapidity spectra may serve as a sensitive test for the hadronic and QGP scenarios 
of strangeness production in p annihilation on nuclei. T he p ossibility for such measurements 
opens up in the planned PANDA experiment at FAIR 10 1. 
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Appendix A: Strangeness production in NN annihilation at rest 

In this Appendix we list the probabihties for the various channels with strangeness for 
pp and pn annihilation at rest (Tables IVl and IVTl respectively). We also provide references 
on source papers. Footnotes explain how some of the probabilities have been obtained from 
experimental data. Many of the partial probabilities for the pp annihilation at rest can be 
also found in the data compilation [41i] . 
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TABLE V: Probabilities of the different channels with strange mesons for pp annihilation at rest 
(in %). 



Channel 


Prob. 


Ref. 


Channel 




Prob. 


Ref. 


K+K- "-^ 


0.091 


[42,43] 


K*+K-^^ 




0.070 


[44J 


^0^0 a) 


0.091 


[42,43j 


K*-K+7r° 




0.070 


m 


K+K*~ 


0.071 


[45j 


^*0^0^0 6) 




0.070 




K-K*+ 


0.071 


m 


^*0^0^0 6) 




0.070 




K^K*^ 


0.060 


[46] 


K*^K~Tr+ 




0.085 


m 


^0^*0 


0.060 


[46] 


K*°K+7r- 




0.085 


m 


K*+K*- 


0.225 


m 


KOk\o^^ 6) 




0.035 


[44,47] 


K*^K*^ 


0.225 


m 


i^+K-^OyrO *>) 




0.035 




K^K^TT^ 1) 


0.146 


[46] 


i^OKOvr+^-vrO 


7) 


0.068 


[48,49] 


K+K-TT^ ^) 


0.146 




i^+K^vr+vr^Tr' 


3 b) 


0.068 




K^K-7r+ 2) 


0.142 


m 


i^OK-^-^+^+ 8) 


0.059 


[48] 


K^K+TT- 2) 


0.142 


m 


^°K+7r+7r"7r" 


- 8) 


0.059 


[48] 


K^K'^ri 3) 


0.050 


[48] 


i^OK-Tr+vrOvrO 


b) 


0.042 




K+K-T] ^) 


0.050 




i^OK+TT-TrOvrO 


b) 


0.042 




K^K^u 3.«) 


0.232 


[50] 


i^+K-TrOvrOvrO 


b) 


0.012 




K+K~Lo '*) 


0.232 


[50] 


i^OKO^OTrOvrO ^) 


0.012 




KOR^pO 4,a) 


0.202 


[44,49,51,52] 


(/)7r+7r^ 




0.054 


[49] 


K+K-p^ 4,a) 


0.202 


[44,49,51,52] 


</>7r07rO ^) 




0.011 




K^K-p+ 5) 


0.234 


[4^ 


<^P° 




0.034 


[49] 


i^OK+p~ 5) 


0.234 


m 


<^7r° 




0.019 


[49] 


iv:*+i?07r- 


0.230 


m 


<^r? 




0.004 


[49] 


K*-ifO^+ 


0.230 


m 


(^o; 




0.030 


[49] 



^^Nonresonance contribution obtained as B^^^o^o = Bk^Ksttq + BklKlttq = 1-46 ■ 10 ^. 
^^Nonresonance contributions obtained as 

rjnonres rjnonres OK/'D D D ^ 

= 0.5 ■ (4.25 - 0.85 - 0.57) ■ 10"^ = 1.42 ■ 10"^ . 
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^^Obtained as Bj^okom = BksKsM + Bk^^KlM-, where M = t], u. 
^^ Obtained as 

BKOROpO = Bj^0ji0^+^- - Bk*±KOtt^,K*±~^K<^tt± - Bk*+K'-,K''+^KO-k+,K*-^K°-k- 



-0.34 ■ {B^^+^- + B^po) = (7.45 - 3.63 - 0.57 - 0.34 • (0.54 + 0.34)) • 10"^ = 2.96 • 10"^ , 

2 

Bx+K-pO = Bk+K-tt+tt- — Bj^*OK*O^K*O^K+'n-,K*O^K-Tr+ ~ -^Bk*Ok±tt^ — 0.49 ■ {B^^^+t,- + B^po] 

= (3.6 - 1.0 - - ■ 1.7 - 0.49 ■ (0.54 + 0.34)) ■ 10"^ = 1.07 • 10"^ . 
^^ Obtained as 

BrOR^p^ = BxOK^n'fTT'' — BK*OK±TT^,K*O^KOnO — B k*^ K±-k^ ,K*^^K°-K'f — Bk*±K°-kT ^K*±^K±-k° 
~Bk*+K'"',K''+^K°tt+{K+tt°),K*-^K-it^{K^^'k-) ~ B x*Ox*0 ^^*0^k°tt°{K+-w-),K*°^K--k+{K°-w°) 
— Bf^,OxOT,OK*O^K'Tr+ ~ Bk''OkOt,Ok*O^K+tt- 

= (10.38 - 0.57 - 0.97 - 1.04 - 0.77 - 1.43 1.4) ■ 10"^ = 4.67 ■ 10"^ . 

^^Nonresonance contribution obtained as 

B'ko'rOt.Ot.o = Bj^OKO^o^o — -{Bk-oRo-kO + Bk*ok»-ko) — -^Bj^,ok,o — Q-MB^t^o^^o 

= (1.1 1.4 2.25 - 0.34 ■ 0.11) ■ 10"^ = 0.35 ■ 10"^ , 



assuming that Bj^oj^o^o^a = B^j^KsnOnO according to |47l |. 
^^Nonresonance contribution obtained as 

-^x°°Fo7r+7r-7rO = ^ k°k07t+7t-7t° ~ 0-28 ■ Bj^o^a^^ — 0.89 ■ Bj^ofio^ — 0.34 ■ (0.28 ■ 5<^^ + 0.89 ■ B^^J) 
= (2.98 - 0.28 • 0.50 - 0.89 • 2.32 - 0.34 ■ (0.28 ■ 0.036 + 0.89 • 0.30)) • 10~^ = 0.68 ■ 10"^ . 

^Obtained as Bj^ok+tt+w-tt- = BxoK-n-TT+TT+ = BksK±it±-k'^itt- 
"^^ Averaged with respect to the charged and neutral K or K* production channels. 
''^Obtained by isospin relations from the branching ratio of another channel. 
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TABLE VI: Probabilities of the different channels with strange mesons for pn annihilation at rest 

(in %). 



Channel 


Prob. 


Ref. 


Channel 


Prob. 


Ref. 


K'^K- 


0.147 


m 


K*OK^TT' ^) 


0.130 


[53j 


KOk*- 1) 


0.067 


M 


K*+K-7r- 8) 


0.154 


[53j 


K-K*^ 1) 


0.067 


M 


K*-K+7r- 8) 


0.154 


[53j 


K*-K*^ 2) 


0.184 


M 


^O^-^O^O a) 


0.043 




i^o^-^o 3) 


0.316 


M 


K^K+TT-TT-TT^ 


0.016 


[53j 


i^OKOvr- 4) 


0.432 


m 


KOK°TT-TT-TT+ 


0.075 


[53j 


K+K-TT- 5) 


0.513 


m 


K-K+7r-7r-7r+ "^ 


0.075 




K^K^T] ") 


0.070 




KOK^TT^TT^TT- ") 


0.052 




K^K-uj 


0.350 


M 


K+K~Tr°TT^TT~ ") 


0.052 




K^K-p^ 6) 


0.150 


[53j 


K^K-tt+tt-tt'' 9) 


0.025 


[53j 


K^K^p- 7) 


0.770 


[53] 


^O^-^O^O^O a) 


0.015 




K+K-p- "^ 


0.770 




C^TT-TT° «) 


0.065 




K*-K^Tr^ 7) 


0.245 


[53j 


cAp-"^) 


0.068 




K*°K~7r° a) 


0.245 




^TT^ 


0.088 


[5^ 


K*^K^TT- 7) 


0.130 


[53] 









^^ Obtained by averaging the transition rates into K*K + K*K from ^5*0 / = 1 and 
^5*1 1 = 1 states as 

Bk*k+k*k = 0.5 ■ (2.0 + 24.6) ■ 10"^ = 13.3 ■ 10"^ , 
Bkok— = Bk-k'O = 0.5 • Bk*k+k*k = 6-7 ■ 10"^ . 

^■'Reconstructed from the partial contribution into the K^K^ti'ti' channel. 
^^Nonresonance contribution calculated as 

Brok-ttO = Bkor-ttO — Bkok*-,K*-^K-tt<^ — Bk-K*o,K''O^KOttO 

= Bkok-^^ - o ■ Bkok'-+k-K'0 = (36.0 - - ■ 13.3) ■ 10"^ = 31.6 ■ 10"^ . 

o o 

^^Nonresonance contribution evaluated as 

nnonres r) _ r) _ r) 

^KOROtt- ~ ^ROROtt- ~ ^ROR—^R — ^ROn- ~ -D(/.7r-,,^^_ft'»Xg 
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= (50.6 - I • 6.7 - 0.34 ■ 8.8) ■ 10"^ = 43.2 ■ 10"^ . 
^^Nonresonance contribution evaluated as 

ononres tj tj tj 

= (60.0 - - ■ 6.7 - 0.49 • 8.8) ■ 10"^ = 51.3 ■ 10"^ . 

^•'Reconstructed from the partial contribution to the K^K^n^n^ channel. 
^•' Reconstructed from the partial contribution to the K^Kgn^n^ channel. 
*•' Reconstructed from the partial contribution to the K'^K^n^n^ and K'^K'^Tc^n^ chan- 
nels. 

^^Nonresonance contribution evaluated as 

nnonrcs rj 75 

-'^KOK-tt+tt-ttO — -DKOR-n+TT-TTO — ^KOK-u},u)^iT+n-iT0 

= (33.6 - 0.89 ■ 35.0) ■ 10'^ = 2.5 ■ 10"^ . 
") Obtained by isospin relations from the branching ratio of another channel. 
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Appendix B: Elementary cross sections 



In this Appendix we describe the cross sections of the hyperon production in A^A^ coUi- 
sions. Other partial cross sections of the NN colhsions are described in Appendix B.2 of 
ref. [ll|. 

The pp — )■ A A and pp — )■ AS° + (c.c.) cross sections (in mb) are parameterized as a 
function of the invariant energy ^/s (in GeV) by the following expressions: 



^pp->AA 



'0, 

0.0357(v^ - 2.232)0-5 + 9.005(v^ - 2.232)l■^ 
^ max(0.01, 0.726(s/4.982 - l)°-^^^(4.982/s)=^-^5°) 



C"pp^ASO+C.C. 



0, 



v/s < 2.232 

2.232 < v^ < 2.252 

2.252 < v^ 

^/s < 2.308 
2.308 < v^ . 



fBll 



(B2) 



0.184(s/5.327- 1)0-437^5 327/^)1.850 ^ 

The functional form in (IBip at lower energies is adopted from ref. 55|, however, with slightly 
different numerical parameters. The parameterization used in fIBip at higher energies and 
in ( ]B2|) is taken from ref. 56| . The numerical parameters in ( IBip and i \B2\i are obtained by 



the fit to the world data on the cross sections pp — > A A and pp — )■ AS° + c.c. from 32 1. 

For the pp — )■ S+S^ cross section only the value of 2 /ib at 3 GeV/c is known experimen- 
tally 32|. Thus, we simply assume the constant cross section above threshold: 

/ib, v^ < 2.630 GeV 
2 /ib, 2.630 GeV < ^ . 

The cross sections involving n, n, A and A in the initial state are obtained from the 
isospin relations as 



a. 



pp^rr.^ 



(^pp^EOEO 



(B3) 



^nn^AA 



'^pp^AA ) 



rB tB. tB rBn tT,\^ 



(^BB^AT.+c.c. = 2cro^/ 4;/ /JI4) 

1=0,1 






_B^ 



(B4) 
(B5) 

(B6) 



where (Tq = cTpp-i. Aso+c.c. 5 ^i = 2crpp_^=+=-, B stands for A^ or A, and only the combinations 
NN , AN and NA of incoming particles are taken into account. Equations (1B5I) and (]B6P 
can be simplified, which gives the following: 



C"nn^ASO+C.C. 



O-Q 



(B7) 
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C^pn^AE-+AS- = 2(To , (B8) 

^A-p^ASO+C.C. = '^AOn^AsO+C.C. = ^0 , (B9) 

C"A-n^AS-+A2- = '^AOp->.AS++AS+ = ^^Tq , (BIO) 

3 

(^A — p-s>AE-+A2- = C"A+n->AS++AE+ = ^(^O (BH) 

'^nn-->S-H++SOsO = (^pn^S-SO = f^l , (B12) 

'7A-p^H-H++H0go = CrAO„_^=-g++H0HO = ^(^l ; (B13) 

C^A-n^H-go = O-A0p^S0S+ = T^"! , (B14) 



^A--p^H-HO = ^A+n^HOH+ = 7^1 (B15) 



3 

— ( 

4 



The partial cross sections for the different outgoing channels shown with a "+" sign in (IB7p - 
(IB 13 II are equal to each other. The expressions for the charged conjugated reaction channels 
are obtained by replacing all particles by the corresponding antiparticles. 

The angular differential cross sections for the NN — t- AA process at ^/s > 2.37 GeV 
{P\sh > 1-830 GeV/c) and for NN — )■ AS, AS processes are chosen according to the phe- 



nomenological Regge-like fit of ref. [57|. At y/s < 2.37 GeV, we fitted the experimental 
angular distributions for the NN — )■ AA scattering of refs. [55|, ISSMGiI] by the following 
expression: 

— ^^ — = [ao + ai exp(a2(tmin - t))] /ib/GeV^ , (B16) 

where tmin = -{Pnn - Paa^ (in GeV^), 

ao = (19.1946 + 132202e^-^)/(l + 42592e^-^) , (B17) 

ai = (277.044e°-2^ - 3145.39e + 4993140e2-^)/(l + 19075.7e^) , (BIS) 

r 67357.8e - 257.692, e < 0.0033 
02 = <^ (B19) 

i 4.66944 ln(1.49962e), 0.0033 < e . 

where e = \^ — 2m\ (in GeV). For the BB — )■ YY processes involving A or A in the initial 
state and for the BB — )■ SS processes we assume the isotropic angular differential cross 
sections in the cm. frame. 
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